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Results are presented from four CDF analyses involving heavy quark production in proton- 
antiproton collisions at center of mass energy 1.96 TeV. The shapes of b-jets are found to be broader 
than inclusive predictions and broader than both PYTHIA and HERWIG defaults. A measurement 
of the production cross section for ip(2S) is consistent with Run 1 results and with theoretical 
predictions associated with parton distribution function energy dependence. The inclusive b-jet 
j production cross section is also consistent with theoretical predictions over six orders of magnitude. 

The bb differential production cross section is compared to several theoretical models and found to 
be best described by MC@NLO + JIMMY. 
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INTRODUCTION 



in 

We report on four measurements by the CDF Experiment [H of processes involving heavy quark production. All 
D ' use data recorded at center of mass energy 1.96 TeV at the Fermilab Tevatron Collider during Run II. 
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II. MEASUREMENT OF fr-JET SHAPES IN INCLUSIVE JET PRODUCTION 

The structure of jets derives from the gluon emissions from the primary parton involved. In heavy quark jets, the 
quark decay must be modeled as well as the underlying event. Multi-gluon emission is difficult to calculate so is 
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often approximated by parton shower models. Jet shape is known to depend upon whether the primary parton is a 
quark or a gluon, and it is also expected to depend upon flavor. It is expected to depend as well upon the production 
mechanism; for example, the b and 6 from gluon splitting are expected to be often in the same jet 0, producing a 
broader jet than does flavor creation. Measurement of jet shapes casts light on all these aspects of jet evolution and 
qq ' heavy flavor production. 

To address the question of whether the fraction of 6-jets originating from gluon splitting, and its evolution, is 
. well described by contemporary models, we define the jet shape, a measure of the fraction of total jet transverse 
momentum px, inside a given radius in the space of rapidity y and azimuthal angle <f>. This is ^(r/R) = ( p^lo^R) ) ' 
the fraction of total px in cone R carried by particles in subcone r, averaged over an ensemble of jets, normalized to 
1, and defined such that ^(0) = and particles outside the cone are excluded. 

The fo-jet event selection begins with triggers based solely on calorimeter transverse energy £x- Trigger Level 1 
requires one calorimetric trigger tower above a threshold between 5 and 10 GeV. Level 2 seeks clusters about the 
Level 1 towers among adjacent towers above 1 GeV and requires that at least one cluster exceed a threshold between 
15 and 90 GeV. Level 3 applies the Run 1 cone algorithm [3| and imposes E™ ln greater than 20 to 100 GeV. The jet 
minimum px requirement is set to ensure that only events for which the trigger is > 99% efficient are used. Offline 
the jets are reconstructed by the Midpoint Cone Algorithm [J]. The selection of b-jets is then enhanced by requiring 
a secondary vertex. Because 6's tend to be found near the jet axis, the search for them uses a cone of radius R = 0.4 
in (y, 4>) space. Tracks are ranked by reconstruction quality, including the contribution from the distance do of closest 
approach to the primary vertex. The algorithm next attempts to reconstruct a secondary vertex, beginning with the 
highest quality track. If this succeeds, a cut is placed on the significance of the two-dimensional projection along 
the jet axis of the distance between the primary and secondary vertices. General event quality cuts follow, including 
a rejection of multiple interactions (there must be one and only one primary vertex with \z\ < 50 cm), a rejection 
of cosmics based on the significance of the missing Et, and a requirement of jet centrality: |j/j C t| < 0.7. The jet 



Pt is corrected to the hadron level by matching hadron and calorimeter level jets in Monte Carlo, a correction that 
increases the pt by 20% to 10% with no change in jet shape. 

The requirement of the secondary vertex tag biases the measured jet shapes by demanding clean, well-defined 
tracks. An r-dependent correction for this is needed. There are additionally corrections for the presence of non-6 jets 
misidentified as 6-jets, corrections for detector effects (for example, the number of calorimeter towers scales with pt), 
and a correction for the presence of a second b quark in the jet. The data are fitted to two-6 and one-6 templates 
separately. The final function to be fitted is given by *£ ad = C had (r/i?) • ^ r / R )-^-P^^-^ R ^^~ b ^/ R \ w here 
C had is the unfolding factor *h ad (r / R)/^ b dct (r / R), is the measured jet shape for the tagged sample, ^d° t n ~ b is 
the measured inclusive jet shape, pb is the purity of b's in the tagged jet sample, and bb and 6 non -b correct for biases 
to the jet shape arising from the secondary vertex requirement. The data are compared to predictions by PYTHIA, 
with Tune A and HERWIG v 6.506 The principle systematics are associated with the choice of Monte Carlo 
model, the effect of the calorimeter model, the jet energy calibration, and assumptions about the c-quark content of 
the jets. 

The results of the study /7j], for integrated luminosity 300 pb -1 , are shown in Figure [IJ for 4 bins inpr< In each bin, 
the data are compared to the PYTHIA prediction for inclusive jets as well as the PYTHIA and HERWIG predictions 
for single b production and for single b production with the default fraction diminished by 20%. One sees that the 
data are incompatible with the inclusive prediction: the jet shapes are evidently influenced by the presence of heavy 
quarks, and &-jets are broader. One also sees that the data jets are broader than the PYTHIA and HERWIG defaults; 
the leading order models underestimate the fraction of 6's from gluon splitting and the fraction of jets with two 6's. 



III. MEASUREMENT OF THE PRODUCTION CROSS SECTION OF THE tp(2S) 

The mechanism for producing heavy vector mesons in hadron collisions is not well understood. CDF Run I data 
on production cross sections for prompt J/ip and ip(2S) were one to two orders of magnitude above predictions by 
color singlet models. Subsequent theoretical efforts 0] with adjusted production matrix elements match the cross 
section but pre dict increasing transverse polarization with production px that is not confirmed by the data Q ■ New 
approaches 1(3, II | have been proposed. The ij)(2S) is a good testing ground for studying charmonium hadropro- 



duction as there are no significant charmonium states above it to produce feeddown. The goal of this study is a 
measurement of a(pp -> i/j(2S)) ■ BR(ip(2S) -> fi + fJ.~) for 2 < p T < 30 GeV/c. 

The data selection begins with muons reconstructed in four layers of the central tracking chamber and matched to 
three to four hit tracks in the muon detector. Events that pass the dimuon trigger, which requires two opposite sign 
tracks, each with significant px, are required to have three hits in the silicon vertex detector SVX II. The ip(2S) mass 
and lifetime are fit with a lifetime function, and requirements are placed on the dimuon mass, 3.5 < m(pfi) < 3.8 
GeV/c 2 , and rapidity \y{p,p)\ < 0.6. 

CDF separates the signal from background, and the prompt from the 6-decay ip's, with an unbinned maximum 
likelihood fit in candidate mass and proper decay leng th ct. The mass separates the signal from the background. The 
signal is modeled with the Crystal Ball Function [l2j , a Gaussian core with a low side tail. The mass background is 
modeled with a first order polynomial. The prompt signal is separated from the feeddown with a ci-fit. The prompt 
signal is given by a double Gaussian centered on zero. The long-lived signal is modeled by an exponential convoluted 
with a Gaussian. The lifetime background is modeled by a prompt double Gaussian plus symmetric, positive-ci, and 
negative-ci long-lived components. The likelihood function is given by L = /s-P s mass (/p-Pp* + (1 — Ip)Pei»g) + (1 — 
fs)PbkT(U^ P sym + f+P+ + f-P C - + (1 - /sym - /+ - f-)P?), where the f's are fractions of signal (s, from the 
total number of candidates in the fit), prompt (p), symmetric long-lived background (sym), positive-rf long-lived 
background (+), and negative-ci long-lived background(-). The P's are probability density functions for signal mass 
(P s mass ), linear mass background (-Pbkg SS ), prompt signal double Gaussian proper time (Pp 1 ), an exponential convolved 
with a Gaussian for &-decay signal (Pe®q), and lifetime backgrounds (-P s C y m , -P+ 4 , and 

The efficiency calculation is given by the product of six terms whose values range from 95.3% to 99.9%. The 




FIG. 1: The measured integrated b-jet shapes for four pt bins. The data are shown as open squares where the error bars 
represent the statistical and total uncertainties. The statistical uncertainties are smaller than the squares. The data are 
compared to PYTHIA Tune A (solid lines) and HERWIG (dashed lines). The uppermost solid curve is the inclusive prediction; 
the next solid curve, and the uppermost dashed curve, are the predictions with the default fn, fraction; and the lowest solid 
curve and lowest dashed curve are the predictions with the default lb fraction diminished by 20%. 



acceptance calculation is complicated by the fact that acceptance depends upon polarization per px bin, yet CDF 
polarization data on the tp(2S) are too weak for direct application. The J/ip polarization, for which strong CDF 
data exist, are, however, expected to correlate highly with prompt ip(2S) polarization, so after confirming that the 
J/ip and ip(2S) polarization data are consistent, CDF applies the J/ip polarization values to the tp(2S) pt bins. 
Monte Carlo i(>(2S) samples are generated with fixed polarizations at the extreme values of and -1 and with flat 
distributions in pr, ?7, and <fi. These samples are simulated in the CDF detector and reconstructed. We compute the 
product of geometrical acceptance and trigger efficiency, A = — — (pt) x (^^(pt)/JV — (pt)) ^ w ] :iere f or Tygen g enera ted 
events, iV rcc survive geometric and reconstruction requirements and N eS survive the trigger. A ranges from 0.0053 
to 0.2548. We interpolate the data to intermediate polarization values predicted from the J/tp data. The principle 
systematic uncertainties on this measurement are associated with the luminosity, the reconstruction efficiency, the 
dimuon trigger efficiency, and the ip(2S) polarization. The variation of the prescale during the course of collecting 
integrated luminosity 1.1 fb _1 produces an effective integrated luminosity of 954 pb _1 . 

The resulting prompt ip(2S) differential cross section is shown in Figure[21 One sees that in all three cases, the Run 
I result is confirmed with an order of magnitude more statistics. The measured value of cr(pp — > tp(2S)) ■ BR(iJj(2S) — > 
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FIG. 2: The prompt production cross section versus pt for data from CDF Run I (triangles) and Run II (circles). 



) is 0.68 ± 0.01 ± 0.06 nb. The ratio of the Run II to the Run I result, for identical ranges in pr, is 18 ± 19%. 
The theoretical prediction [131 ] for this ratio, including the energy dependence of the parton distribution functions, 
is 14 ± 8%. 



IV. STUDY OF INCLUSIVE b-JET PRODUCTION 

Measurements of b quark production in hadronic collisions test perturbative QCD. The CDF Run I inclusive 
B meson cross section measurement [l4| motivated theoretical developments beyond next-to-leading order 
Continued comparisons, at higher px and with increased statistics, may motivate further improvements. 

This study considers events whose primary vertex has \z\ < 50 cm, whose missing is not significant, and whose 
kinematical variables lie in the range 38 < p£ < 400 GeV/c and |y jct | < 0.7, with jets reconstructed with the 
Midpoint Cone Algorithm and heavy flavor jets tagged via their secondary vertices. The jet energy corrections are 
based on minimum bias events which measure the pt deposited in the calorimeter, as a function of the number of 
primary vertices, and indicate a value of 0.93 ± 0.14 GeV/c per extra vertex. Monte Carlo events from PYTHIA 
6.203, with Tune A and the CTEQ5L parton distribution function set, are used to correct the measured energy for 
effects of partially instrumented regions and calorimeter nonlinearities. Typical unfolding factors range from 1.6 to 
2.1. The b quark content of the jet is inferred from the shape of the invariant mass of all charged tracks attached to 
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FIG. 3: The measured secondary vertex mass distribution for jets with corrected jet pt in the range - 
to the result from the fit to Monte Carlo templates for b and non-6 jets. 



to 92 GeV/c, compared 



the secondary vertex, as shown in Figure [3] The measurement is compared to the theoretical prediction described in 
0. 

The measured differential cross section, d f b ~J at = N tass <.dh j g s j lown m Figure g] for integrated lumi- 

aPTO-y eb-tagAj/J * Ap J T J Ldt 

nosity 300 pb . Here -/Vt agg ed is the number of tagged jets in eachpr bin, Apx is the bin size, fb is the fraction of jets 
in the tagged sample (this ranges from 0.35 at low pt to 0.14 at high), £b-tag is the 6-tagging efficiency, Ayi et is the 
jet rapidity range, and L is the luminosity. The measured cross section is consistent with the theoretical prediction 
throughout the full six orders of magnitude of range in p^. It exhibits a strong dependence on scale /i, however, 
which suggests that higher orders may make large contributions. The largest systematic uncertainties derive from 
the luminosity, the jet energy scale, the jet energy resolution, the unfolding process, the 6-tagging efficiency, and the 
fraction of 6-jets, and produce a total systematic uncertainty that ranges from ±25% at low pp to ±70% at high pr- 



V. STUDY OF THE PRODUCTION OF bb DUETS 

Momentum conservation requires that the azimuthal angle <fi between a b and b produced at lowest order must 
be 180°. Higher order QCD processes produce additional partons in the final state, modifying the range of allowed 
azimuthal angle difference, A(f>. Order diagrams are thought to contribute the same magnitude to the cross section 
as order a 2 s diagrams. Measuring the cross section as a function of Acf> provides information on the contributions of 
the leading order and next-to-leading order terms. 

The event selection requires two central calorimeter towers with significant Et in association with two tracks in the 
central tracking chamber with significant px- The calorimeter towers arc clustered, and the tracks are reconstructed 
with the Silicon Vertex Trigger Events having two clusters with significant £r associated with two tracks 

having significant displaced vertices are retained. Jets are reconstructed with the Run I cone algorithm. One or 
more high-quality primary vertices are required. Corrections are applied to the jet energy scale as well as to the 
energies of 6-jets (whose fragmentation, for example, is harder than that of light quarks). The presence of two jets 
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FIG. 4: The measured inclusive 6-jet cross section (filled circles) as a function of p 1 ^ compared with the next-to-leading order 
perturbative QCD prediction (empty crosses). The shaded band represents the systematic uncertainty on the data and the 
dashed band represents the uncertainty on the theory. 



with significant E^, both 6-tagged and with displaced secondary vertices, is confirmed. The data are compared to 
predictions by PYTHIA, with Tune A; HERWIG; and MC@NLO 18] with the HERWIG parton shower and the 
underlying event generated by Jimmy 4.3 [l9j |. The systematics are dominated by the jet energy scale and range 
from 20 to 30%. 

For integrated luminosity 260 pb -1 , all three Monte Carlos show similar agreement with data distibutions of bb 
production cross section versus leading jet Et and versus dijet invariant mass. The cross section as a function of A</>, 
however, is significantly better modeled by MC@NLO, as shown in Figure [51 The peak at large angles reflects flavor 
creation, and the excess at small angles is due to higher order diagrams and multiple interactions. The measured 
total cross section for |r? 12 | < 1.2, E T1 > 35 GeV, and E T2 > 32 GeV is a = 5664 ± 168 ± 1270 pb. 
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FIG. 5: The bb jet cross section as a function of the dijet A<j) correlation. Data are compared to hadron level cross sections 
obtained using MC@NLO+JIMMY, PYTHIA, and HERWIG+ JIMMY. The shaded area represents the total systematic 
uncertainty on the data. 
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